Quantitative comparison of single- and two-particle properties in the cuprates 
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We explore the strong variations of the electronic properties of copper-oxygen compounds across 
the doping phase diagram in a quantitative way. To this end we calculate the electronic Raman 
response on the basis of results from angle-resolved photoemission spectroscopy (ARPES). In the 
limits of our approximations we find agreement on the overdoped side and pronounced discrepancies 
at lower doping. In contrast to the successful approach for the transport properties at low ener- 
gies, the Raman and the ARPES data cannot be reconciled by adding angle-dependent momentum 
scattering. We discuss possible routes towards an explanation of the suppression of spectral weight 
close to the (7r, 0) points which sets in abruptly close to 21% doping. 

PACS numbers: 



I. INTRODUCTION 

Superconductivity in the copper-oxygen compounds 
(cuprates) 1 occurs at transition temperatures T c above 
100 K with the origin still controversial. On a basic level 
the problem can be cast into the question as to whether 
the condensed state at low-temperature is an instabil- 
ity of the normal metal as in conventional systems or a 
new ground state. If superconductivity results from in- 
teractions between the electrons via an intermediary bo- 
son the theoretical approach proposed by Eliashberg 2 4 
should provide an at least qualitative description of the 
condensation of electrons into Cooper pairs. Signatures 
of the coupling boson should also determine the prop- 
erties above T c . Traditionally, one looks at renormaliza- 
tion effects in the single-particle or transport responses as 
a function of temperature, specifically around T c where 
abrupt changes are to be expected 5 7 . If the spectral 
shape of the coupling function can be retrieved 8 10 one 
may be in a position to derive the relevant interaction(s). 
However, the results in the cuprates are far from con- 
verging into a unified picture. The unsuccessful search 
for clear signatures of a dominating retarded interaction 
was one reason why new ground states were and are being 
studied intensively 11 13 . 

There were early proposals that the superconductiv- 
ity in the cuprates could be a property of the electrons 
alone. A prominent example is the resonating valence 
bond (RVB) state with antiferromagnetic coupling of 
nearest neighbor spins 11 which emerges from the Neel 
state upon doping 12 . From an experimental point of view 
this proposal is much harder to pin down than supercon- 
ductivity as an instability or a small perturbation of the 
normal state. There is indeed an arsenal of methods to 
analyze data in terms of an Eliashberg-type strong cou- 
pling approach 3 ' 7 10 ' 14 while ideas for experiments un- 
veiling new ground states are scarce. 



Our Raman experiments aim at disentangling strong 
coupling effects evolving mainly as a function of 
temperature 9 from phase transitions and cross-over phe- 
nomena occurring as a function of doping (see also 
Ref. 15 ). We develop an analysis using a phenomenolog- 
ical description of the response 16 . Starting from single- 
particle properties observed by angle-resolved photoemis- 
sion spectroscopy (ARPES) we calculate the electronic 
Raman response (ERS) in both the normal and the su- 
perconducting states. While ARPES measures occupied 
single-electron states ERS yields - similarly as, e.g., op- 
tical transport (IR) - a weighted convolution of occu- 
pied and empty states. Therefore the two-particle spec- 
tra contain additional information originating from the 
interaction between the hole and the electron created in 
the scattering process. 



II. EXPERIMENT AND THEORETICAL 
BACKGROUND 

In the Raman process an electron is scattered through 
an intermediate high-energy state into an unoccupied 
level (for a more detailed discussion see Refs. 15 ' 17 ). 
The energy transferred to the electron corresponds to 
the difference between the incoming and the outgoing 
photon, Q = uji — uj s . As a result of the Coulomb 
interaction only charge fluctuations inside the unit cell 
can be observed. The resulting selection rules lead to 
form factors 7^ weighing out symmetry-specific regions 
of the Brillouin zone. For crossed photon polarizations 
aligned along the principle axes (of the Cu02 plane) the 
B2g symmetry is projected out and r )^ 2g cx sin k x sin k y . 
In Bi g symmetry with both polarizations rotated by 45° 
the form factor reads 7jf lp cx cos k x — cos k y as shown on 
the r.h.s. of Fig. 1. Excitations with A\ g symmetry will 
not be considered here. 
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For the B symmetries the Raman response function 
x" 7 (fi) can be expressed directly in terms of a generalized 
susceptibility Xa,b(fy which, in Nambu representation, 
reads 17 
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with T the temperature, TV the number of k-points a, b 
the vertices and G(k, iuj) the matrix Green function. The 
vertices can be expressed in terms of Pauli matrices f^. 
For charge excitations, a(k) = f^a{k). a(k) can be either 
unity or one of the vertices 7^ or the electron velocity 
to yield the density or the Raman response or, respec- 
tively, the conductivity. In general, the vertex b must 
be renormalized 17 but in the approximation here we use 
only bare vertices. G'/^k, uj) = — 7rA(k, uj) is the imag- 
inary part of the renormalized Green function which is 
proportional to the spectral function A(k, uj) which, for 
the occupied states, can be measured by ARPES. 

In strong-coupling theory the components of G(k, uj) 
read 7 ' 19 
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etc. with = — /i the bare band structure. The 
complex functions Z(k,cj), x(k,ct;), and <£(k, co>) must be 
found self consistently and describe all interactions of the 
electrons. In the weak coupling limit Z = 1, x = 0? and 
$ = Ak- For Z / 1 and x finite the usual self energy 
E = £' + z£" is related to Z and x as E" = -cjZ" + 
x //7,20,2i 5 z = 1 - E/cj and Z = (1 + E'/cj^below 
and above T c , respectively, and $ = ZAk 20 in lowest 
order approximation. Since E' is ant i- symmetric and E" 
vanishes faster than uj below T c even in the presence of 
impurities or nodes, Z is always defined. In the normal 
state vanishes and G± \ = Gq 1 — E. 

For the calculation of x" 7 (^) we need both the oc- 
cupied and the empty states and model functions for 
A(k, uj) must be derived from the experimental ARPES 
spectra A exp (k, uj) being cut off at the chemical poten- 
tial /i(T) by the Fermi function f(T,uj). Ideally, one has 
to microscopically derive expressions for the bare band 
structure £k = <^c — [i and for E, Z, x, and <E> in order to 
obtain a spectral function A(k, uj) which, after multipli- 
cation with the matrix element |M/^| 2 and convolution 
with the experimental resolution function R exp , repro- 
duces A exp = (/ A\M f ,i\ 2 ) <g> # exp . 

Since we do not focus on microscopic models here nei- 
ther for the normal nor for the superconducting state 



we use phenomenological approximations for both the 
self energy and the band structure which reproduce the 
observed ARPES data satisfactorily. A tight binding de- 
scription with nearest and next-nearest neighbor hopping 

£k = — 2t(cos k x + cos k y ) + At' (cos k x cos k y ) — fi (4) 

using t = 0.25 eV for the nearest neighbor hopping inte- 
gral and t'/t = 0.35 for La 2 _ x Sr x Cu0 4 (LSCO), t' /t = 
0.40 for Bi 2 Sr 2 CaCu 2 8+ 5 (Bi-2212), and a self energy 
of the form 



E" = - 



+ eg + c k 



(5) 



which is inspired by (but not directly drivable from) 
the marginal Fermi-liquid (mFL) approach, E" = 
— max(a|co>|, /3T) 22 , yield good agreement with photoe- 
mission. The chemical potential fi(T) is adjusted to give 
the correct filling and a momentum dependent elastic 
scattering term Ck is additionally included. For Co = 0, 
E" is strictly linear at T = 0. With Co finite E" varies 
always quadratically at low energy. We assume that E" 
saturates above an energy ujq which approximately corre- 
sponds to the Joffe-Regel limit equivalent to a mean free 
path £ « a with a the lattice constant. The real part of 
the self energy E' is obtained from Eq. (5) by Kramers- 
Kronig transformation. In the superconducting state we 
follow the phenomenology of Inosov et ai 23 . 



III. RESULTS AND DISCUSSION 

We have measured electronic Raman spectra of Bi- 
2212 and LSCO above and below the "critical" doping 
p = x^0.21at which the B\ g response changes quali- 
tatively. In Bi-2212 we have evidence that the transition 
is quite abrupt 24 ' 25 . Neither in Tl 2 Ba 2 CaCu0 6 +5 (Tl- 
2201) nor in LSCO the doping levels are sufficiently dense 
to allow this conclusion. However, the crossover in the 
range of 21% is clearly seen and seems to be one of the 
generic features of the Raman response in the cuprates 15 . 
We now try to answer the question in which way this 
feature can be related to other properties of the cuprates 
as observed for example in ARPES or transport experi- 
ments. 

We start with non-superconducting LSCO (x = 0.26) 
and parameterize the self-energy as a = 1.1, j3 = 2.5, 
co/t = 0.04, c k = 0, and fi(100K)/t = 1.3 to obtain 
the proper filling and satisfactory agreement with the 
ARPES spectra 26 28 . In addition, the temperature de- 
pendence of the resistivity is well reproduced. With these 
parameters, optimized for Lai.74Sro.26Cu04, and using 
Eq. (1) we arrive at the Raman spectra shown in Fig. 1. 
The agreement is remarkable since only a single intensity 
is set at 150 K. Symmetry and temperature dependence 
follow then within the model. 

It has been shown before that the longitudinal and 
the Hall resistivity can be reproduced in a Boltzmann 
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FIG. 1: Temperature dependence of the electronic Raman response of overdoped non-superconducting Lai.74Sro.26Cu04. 
The smooth lines represent the phenomenology described in the text and are based on an analytic approximation to the 
photoemission results. 



approach with similar assumptions 29 . Here, we use full 
k sums to calculate the spectra at finite energies. The 
limitation to the Fermi surface would suppress essential 
structures of the spectra originating from the anisotropy 
of £k and the proximity of the van Hove singularity to ja. 
The agreement reaches beyond the energy range shown 
here and leaves only little spectral weight unexplained up 
to 1 eV. The agreement between various single- and two- 
particle probes at high doping demonstrates that the self 
energy alone captures the essential many-body physics 
and that the lowest order approximation given in Eq. (1) 
is sufficient for the calculation of response functions. 

The only type of microscopic models for E which quali- 
tatively reproduce the variation with energy and temper- 
ature compatible with Eq. (5) are based on fluctuations. 
They lead to "marginal" behavior with strongly reduced 
quasiparticle weight, Z — >• 0, at the Fermi surface, in 
certain doping and temperature ranges. The marginal 
Fermi-liquid model (mFL) 22 was the first proposal which 
predicted expressions similar to Eq. (5). Later, circu- 
lating orbital currents were proposed to be the micro- 
scopic origin of this phenomenology 30 which may leave 
an imprint on the Bragg peaks in neutron scattering 
experiments 31 . Fluctuations of an incipient charge den- 
sity wave (CDW) 32,33 or a Fermi surface deformation 34 
are also candidates. They lead to Fermi liquid like varia- 
tions of E at low temperature, actually well below T c , and 
predict a strong momentum dependence of the self energy 
close to the quantum critical point where the transition 
temperature to the ordered or partially ordered state ap- 
proaches or extrapolates to zero. In fact, indications of 
CDW fluctuations, nematic or even long-ranged ordering 
have been observed in various experiments 35 40 . 

At high doping, p > 0.21, an anisotropy of the self en- 
ergy is not observed. Upon decreasing the doping level 
the ARPES spectra show only little variation and even 
around optimal doping there are well defined quasipar- 



ticle peaks on the entire Fermi surface. All changes ap- 
pear to be continuous 26 28 . In contrast, the B\ g Raman 
spectra change abruptly close to p = 0.21 24 ' 25 . Since 
there is no discontinuity in the ARPES spectra around 
p = 0.21 we can use the same model for E and £k on 
either side of the crossover. The comparison between the 
prediction on the basis of the ARPES results and the 
observed B\ g Raman spectra for samples above and be- 
low p = 0.21 are shown in Fig. 2. In I>2g symmetry, the 
spectra are well described at both doping levels and, be- 
yond that, down to the lowest carrier concentration inside 
the superconducting dome 15 . The B\ g spectra, however, 
drop considerably below the simulation in the range be- 
low 1000 cm -1 . While the spectral shapes are similar in 
Bi-2212 and LSCO the overall intensity in LSCO is sub- 
ject to variations of the cross section due to resonance 
effects. 

The main result here is twofold, (i) A strong inter- 
action sets in below approximately 21% doping which 
(ii) manifests itself predominantly in the two parti- 
cle properties. Indications of this interaction have 
been seen a long time ago by nuclear magnetic reso- 
nance (NMR) 41 which measures a similar response func- 
tion [see Eq. (1)]. With Raman scattering we can 
additionally pin down the momentum dependence 24 . 
In IR experiments anomalies are mainly observed in 
the c-axis conductivity of i?Ba2Cu306+ x (#-123 with 
R = Y, Nd, La) compounds 42 44 while the transport in 
the a — b plane is essentially captured by a Fermi liq- 
uid or mFL phenomenology 45 47 . This is not unexpected 
and compatible with the Raman results when the BZ 
projections of the respective methods are taken into ac- 
count. Then, due to the band structure of the cuprates 48 
in-plane IR and B^ g Raman project predominantly the 
nodes while out-of-plane transport and B\ g Raman are 
more sensitive at the anti-node 49 . 

It has been suggested that an increasing contribu- 
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Raman shift Q (cm 1 ) 

FIG. 2: Doping dependence of the electronic Raman response of E^S^CaC^Os+s and La2- x Sr x Cu04. There is little variation 
for the Big spectra. The B\ g responses change abruptly at p c± 0.21 and cannot be described any more by spectral functions 
derived from the ARPES results in lowest order approximation [Eq. (1)]. 



tion from anisotropic elastic scattering of electrons leads 
to the dichotomy between nodal (7r/2,7r/2) and anti- 
nodal (7r, 0) particles 47 ' 50 ' 51 . We check this proposal 
here by adding an energy and temperature independent 
anisotropic scattering rate varying as 

Ck = ci + Y"( cos k x — cos ky) 2 . (6) 

With c\ and C2 properly selected the variation of the scat- 
tering rate on the Fermi surface as used for Tl-2201 46 ' 52 
(but not any more for LSCO 47 ) is reproduced to within 
a few percent. In Fig. 3 we compare the resulting Raman 
spectra with the experimental data of LSCO at optimal 
doping. For our choice of parameters the phenomenolog- 
ical curves fit the shape and the temperature dependence 
of the Big data reasonably well but deviate now signif- 
icantly from those in B^ g symmetry. In LSCO a fit to 
the Big results is still possible since the temperature de- 
pendence of the spectra remains metallic. In contrast, 
the spectra in Bi-2212 exhibit even a slightly insulating 
variation with temperature around optimal doping 24 and 
agreement between the calculated and the experimental 
Big spectra cannot be obtained any further. We conclude 
that an elastic term is insufficient to reconcile single- and 
two-particle properties at optimal doping in the same 
approximation as at high doping and that dynamic in- 
teractions are at the origin of the observed effects. In 
addition, it is not completely surprising that higher or- 
der corrections may become necessary when entering the 
low doping range. 

The reason why the transport data can be repro- 
duced reasonably well in terms of an isotropic Boltzmann 
approximation 47 is not immediately obvious. The linear 
dispersion, ek ~ ^k(p — Pf), around the Fermi momen- 



tum pf and the specific variation of the current vertices 
used in that approach may be part of the explanation. In 
earlier studies of the superconducting response we actu- 
ally found substantial differences between Fermi surface 
and k integration. After all, the transport is mainly due 
to nodal quasiparticles which survive down to at least 
p = 0.05 and may be mostly blind for the effects we find 
around (-zr, 0). 

The most remarkable results here are the continuous 
versus abrupt changes in ARPES and transport studies 
and, respectively, in the Raman spectra across the critical 
doping level of p « 0.21. Apparently, we encounter a 
transition from an essentially conventional metallic state 
to one of strongly interacting electrons. The origin of 
the dichotomy between single- vs two-particle properties 
has to remain open at the moment. From the Raman 
scattering point of view a doping-dependent elastic term 
can be excluded as an explanation. 

Below T c an additional complication arises. While 
the B2 g response [Fig. 4 (a)] is universal, the Bi g 
spectra clearly reflect sample specific behavior and the 
superconductivity-induced features exhibit a statistically 
significant sample dependence at a given doping level as 
shown in Fig. 4 (b). More quantitatively, the B^ g re- 
sponse [Fig. 4 (a)] scales with the individual T c , and is 
satisfactorily described already by the weak-coupling pre- 
diction. There are only small differences in the intensity 
close to the peak maximum which can be traced back to 
variations in the impurity concentration 59 . In contrast, 
the Big spectra [Fig. 4 (b)] clearly reflect sample spe- 
cific behavior and, beyond that, do not scale with T c but 
rather as (1 — p) 17 ' 53 58 > 60 > 61 . Moreover, on the basis of 
results obtained with applied pressure 62 we conclude that 
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Raman shift Q (cm" 1 ) 

FIG. 3: Temperature dependence of the electronic Raman response of optimally doped Lai.ssSro.isCuCU. The smooth lines 
represent the calculations. Here, an additional anisotropic contribution from elastic collisions according to Eq. (6) with Co = c\ — 
and C2 = O.lt. The inset in (a) shows the spectral function at the Fermi energy [A(k, u — 0)] in false-color representation (red 
high, blue low intensity) for these parameters. The suppression of intensity around (ir, 0) is much stronger than experimentally 
observed at this doping level as shown in panel (c) (From Ref. 27 with permission.) 




B c 

FIG. 4: Electronic Raman response of Bi 2 Sr2CaCu208+5 well below T c close to optimal doping. The energy is given in units 
of ksT c . The hatched area indicates the energy range, where superconducting and normal state spectra merge. The dashed 
and the full smooth (red) line represents the weak and strong coupling predictions, respectively, for the superconducting gap 
structures. 
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internal strain induced by quenched disorder 63 leads to 
the variation of the B\ g spectra below T c . This experi- 
mental fact escaped attention so far but may shed light 
on the origin of the B\ g spectra. It is indeed hard to 
understand their origin in terms of a pure pair-breaking 
effect. 

This qualitative reasoning is fully corroborated by the 
simulations derived from the ARPES spectra. Here, 
Eq. (1) must be formulated for superconductivity, and 
the diagonal (normal) and the off-diagonal (anomalous) 
parts of the Green function 17 are needed. The self ener- 
gies have to be derived in the spirit of a strong coupling 
approach 7 ' 9 ' 20 ' 23 . The agreement with the B 2g spectra 
is comparable to that of the weak-coupling prediction 
(Fig. 4 (a)). However, there is no agreement between the 
simulations and the B\ g Raman spectra [Fig. 4 (b)] al- 
though sample Bi-UD92 came from the same source as 
the one used by Inosov et al. 23 for ARPES. Results for 
samples from other sources lie more or less inside the 
pair-breaking peak derived from ARPES but, similarly 
as in the normal state, the overall experimental intensity 
is generally too small. The larger intensity may originate 
in the Fermi liquid like variation of the self energy pro- 
posed by Inosov et al. 23 to fit the ARPES spectra which 
leads also to discrepancies in B 2g symmetry above the 
maximum [Fig. 4 (a)]. At lower energies E" appears to 
vary closer to linearly 47 ' 59 rather than quadratically. 

IV. CONCLUSIONS 

The quantitative comparison of ARPES and Raman 
spectra in lowest order approximation yields quantitative 
agreement of the normal-state response above p ~ 0.21. 
The different spectral shapes of the B\ g and B 2g sym- 
metry are obtained with a momentum independent self 
energy and can be traced back to the band structure 
alone, specifically to the proximity of the van Hove sin- 
gularity to the chemical potential. At p < 0.21, discrep- 
ancies between the single- and two-particle responses are 
observed predominantly in B\ g symmetry while nodal 
ARPES and B 2g Raman spectra remain consistent by 
and large. The B\ g Raman spectra in both the nor- 
mal state and the superconductivity-induced features are 
progressively suppressed in a frequency range of at least 



1500 cm -1 . As opposed to transport measurements 29 ' 47 , 
a momentum dependent constant relaxation term alone 
which, for instance, could originate from correlated scat- 
tering centers 50 is insufficient to explain the differences. 
We rather conclude, that the origin of the discrepancies 
is due to dynamic processes. The abrupt onset of renor- 
malization effects in the B\ g spectra suggests a crossover 
controlled by a quantum critical point close to the center 
of the superconducting dome 24 . A new type of interac- 
tion appears which is not fully described by the single- 
particle self energy. It is an interesting question as to 
whether or not the suppression of spectral weight has its 
origin in the the same interactions as those which favor 
a new ground state and induce superconductivity. 

Similarly as in the normal state, the B 2g responses be- 
low T c look rather conventional and scale with the indi- 
vidual transition temperatures in the entire doping range 
studied 61 . The explanation in terms of a weak-coupling 
approach is limited to Q < 2Ao. In a strong coupling 
approximation with the self energy and the band struc- 
ture derived from the ARPES results an energy range 
of at least a few A can be assessed. This is not the 
case in B\ g symmetry. At optimal doping, p = 0.16, 
the pair-breaking features observed by Raman scattering 
reflect properties of the individual samples and are in- 
consistent with the single-particle results. This implies 
that the B\ g spectra do not directly and exclusively re- 
flect the maximal gap A . Rather, pairing correlations 
induce features close to 2Ao and, additionally, seem to 
activate another excitation which is particularly clearly 
seen for sample Bi-OPT94. In general, the B\ g response 
appears to provide crucial yet unexplained information 
for the understanding of the relevant interactions in the 
cuprates. 

Acknowledgements 

We are indebted to L. Benfatto, S. Caprara, T.P. Dev- 
ereaux, C. Di Castro, M. Grilli, B. Moritz, and A. Zawad- 
owski for valuable discussions. The work has been sup- 
ported by the DFG via Research Unit FOR 538 (grant- 
no. Ha2071/3). B.M. and R.H. gratefully acknowledge 
support by the Bavarian Californian Technology Center 
(BaCaTeC). 



Permanent address: Mettler- Toledo (Schweiz) GmbH, 

8606 Greifensee, Switzerland 
t Electronic address: hackl@wmi.badw.de 
* Present address: MPI fiir Festkorperforschung, Heisen- 

bergstr. 1, 70569 Stuttgart, Germany 

1 J.G. Bednorz and K.A. Miiller, Z. Phys. B 64, (1986) 189. 

2 G. M. Eliashberg, Zh. Eksperimen. i Teor. Fiz. 38, (1960) 966 
[Sov. Phys. JETP 11, (1960) 696]. 

3 J. P. Carbotte, Rev. Mod. Phys., 62, (1990) 1062. 

4 J. P. Carbotte and F. Marsiglio, The Physics of Superconductors 
(Springer, Berlin 2003) 233. 



5 S. Engelsberg and J.R. Schrieffer, Phys. Rev., 131, (1963) 993. 

6 P. B. Allen and M. L. Cohen, Phys. Rev. Lett., 29, (1972) 1593. 

7 T. Cuk et a/., phys. stat. sol. (b) 242, (2005) 11. 

8 W. L. McMillan and J. M. Rowell, Phys. Rev. Lett. 14, (1965) 
108. 

9 J. Carbotte, E. Schachinger & D. N. Basov, Nature 401, (1999) 
354. 

10 A. A. Kordyuk et al, Eur. Phys. J. ST (2010) in print 

11 P. W. Anderson, Science, 237, (1987) 1196. 

12 Patrick A. Lee, Naoto Nagaosa, Xiao-Gan Wen, Rev. Mod. 
Phys. 78, (2006) 17. 



7 



13 P. W. Anderson, Science, 317, (2007) 1705. 

14 T. Dahm et al, Nature Physics, 5, (2009) 217. 

15 B. Muschler et al, Eur. Phys. J. ST (2010) in print; preprint at 
arXiv:1008.5316 

16 F. Venturini, Doctoral Thesis, Technische Universitat Miinchen 
(2003). 

17 T. P. Devereaux and R. Hackl, Rev. Mod. Phys. 79, (2007) 175. 

18 D. J. Scalapino, Superconductivity I edited by R. D. Parks (Mar- 
cel Dekker, New York 1969) 449 

19 G. Mahan, Many-Particle Physics, (Kluwer Academic/Plenum 
Publishers, New York 2000). 

20 M. Eschrig and M. R. Norman, Phys. Rev. B 67, (2003) 144503. 

21 A. V. Chubukov and M. R. Norman, Phys. Rev. B 70, (2004) 
174505. 

22 C. M. Varma et al, Phys. Rev. Lett. 63, (1989) 1996. 

23 D. Inosov et al, Phys. Rev. B 75, (2007) 172505. 

24 F. Venturini et al, Phys. Rev. Lett. 89, (2002) 107003. 

25 S. Blanc et al, Phys. Rev. B 80, (2009) 140502(R). 

26 A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 
75, (2003) 473. 

27 T. Yoshida et al, Phys. Rev. B 74, (2006) 224510. 

28 J. Chang et al, Phys. Rev. B 78, (2008) 205103. 

29 S. Nakamae et al, Phys. Rev. B 68, (2003) 100502. 

30 CM. Varma, Phys. Rev. B 55, (1997) 14554. 

31 Y. Li et al, Nature 455, (2008) 372. 

32 S. Caprara et al, Phys. Rev. B 59, (1999) 14980. 

33 M. Grilli et al, Physica B 404, (2009) 3070. 

34 W. Metzner, D. Rohe, and S. Andergasen, Phys. Rev. Lett. 91, 
(2003) 066402. 

35 J. M. Tranquada et al, Nature 375, (1995) 561. 

36 J. M. Tranquada et al, Nature 429, (2004) 535. 



37 H.-H. Klauss et al, Phys. Rev. Lett. 85, (2000) 4590. 

38 J. Fink et al, Phys. Rev. B 79, (2009) 100502. 

39 V. Hinkov et al, Eur. Phys. J. ST (2010) in print. 

40 M. Vojta, Eur. Phys. J. ST (2010) in print. 

41 H. Alloul, T. Ohno, and P. Mendels, Phys. Rev. Lett. 63, (1989) 
1700. 

42 C. C. Holmes et al, Phys. Rev. Lett. 71, (1993) 1645. 

43 C. Bernhard et al, Phys. Rev. Lett. 80, (1998) 1767. 

44 D. N. Basov and T. Timusk, Rev. Mod. Phys. 77, (2003) 721. 

45 D. van der Marel et al, Nature 425, (2003) 271. 

46 N. E. Hussey, Eur. Phys. J. B 31, (2003) 495. 

47 R. A. Cooper et al, Science 323, (2009) 603. 

48 O. K. Andersen et al, J. Phys. Chem. Solids 56, (2009) 1573. 

49 T. P. Devereaux, Phys. Rev. B 68, (2003) 094503. 

50 A. P. Kampf and T. P. Devereaux, Phys. Rev. B 56, (1997) 
2360. 

51 E. Abrahams and C. M. Varma, PNAS 97, (2000) 5714. 

52 M. Abdel-Jawad et al, Phys. Rev. Lett. 99, (2007) 107002. 

53 C. Kendziora and A. Rosenberg, Phys. Rev. B 52, (1995) R9867. 

54 R. Nemetschek et al, Phys. Rev. Lett. 78, (1997) 4837. 

55 S. Sugai and T. Hosokawa, Phys. Rev. Lett. 85, (2000) 1112. 

56 M. Opel et al, Phys. Rev. B 61, (2000) 9752 . 

57 F. Venturini et al, J. Phys. Chem. Solids 63, (2002) 2345 . 

58 R. Hackl et al, Adv. in Solid State Physics 45, (2005) 227 . 

59 T. P. Devereaux, Phys. Rev. Lett 79, (1995) 175. 

60 M. Le Tacon et al, Nature Physics 2, (2006) 537 . 

61 N. Munnikes et al, preprint at arXiv:0901.3448. 

62 A. F. Goncharov and V. V. Struzhkin, J. Raman Spectrosc. 34, 
(2003) 532. 

63 H. Eisaki et al, Phys. Rev. B 69, (2004) 064512. 



